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In a recent experimental study Q of the reaction 
^^F(p, 2p)^®0 two-proton events were measured from excita¬ 
tions near a 1“, -E* = 6.15 MeV state in ^®Ne. We calculate 
by means of i?-matrix theory the resonant two-proton pro¬ 
duction cross section and branching ratios. We conclude that 
it is unlikely that two-proton production via population of 
the 1“ state is sufficient to explain the observed two-proton 
events. Alternative sources of such events are discussed. 


of Goldansky that protons should evenly share the energy was 
conhrmed in these studies. It was shown in ^ that the en¬ 
ergy distributions in democratic decays can be described by 
the expression 

dN/de oc Ve(E 2 p - e) |A|2 (1) 

where e is relative energy of the protons, E 2 p is the energy of 
the resonance relative to the 2p threshold, and amplitude A 
depends only weakly on e. 


Introduction. In recent paper Q, states in ^®Ne were popu¬ 
lated in the ^^F-|-p reaction. High quality results for the one- 
proton excitation function are obtained in the energy range 
0.4-2.45 MeV, and the question arises whether exotic two- 
proton emission processes can be seen in such an experiment. 
Because there are no intermediate states in ^^F available for 
sequential decay, the region of the ^®Ne spectrum below 6.5 
MeV might be envisaged as a good place to study simulta¬ 
neous two proton decay. This needs a theory of two proton 
decays. 

We therefore review the theoretical foundations of the two 
mechanisms (as mentioned in the introduction of paper 0): 
(i) diproton mechanism Q and (ii) democratic decay In 
his often-referenced paper Q, Goldansky pointed that a cu¬ 
rious quantum mechanical effect is possible in some proton- 
rich nuclei: what he calls “true two-proton emission” takes 
place if there are no two-body decay channels at all due to 
energy conditions in the subsystems. The Hngerprint of this 
effect should be “the energy correlation between the protons 
during the two-proton decay, which leads to their energies 
being almost equal”. Goldansky also noticed that the esti¬ 
mated penetrability for two s-wave protons is very close to 
the penetrability for the ‘diproton’ (charge 2 ‘particle’ with 
zero energy of internal motion). We remark that this simi¬ 
larity of penetrabilities has led to the widespread idea that 
Goldansky expected two-proton decay to produce two pro¬ 
tons with almost coincident velocities, whereas he predicted 
merely coincident energies. 

A slightly different view on the phenomenon is provided by 
the concept of democratic decay [H-^: it was shown experi¬ 
mentally for decay of the ®Be g.s. that the energy distribution 
between the protons is broad (‘democratic’) and no ‘diproton’- 
type correlation was observed. It was suggested that the rea¬ 
son for the broad distributions is the absence of narrow states 
in all subsystems in this decay. This is less stringent condition 
than the condition for ‘true two-proton decay’ in the paper of 
Goldansky, namely that the “positive binding energy of the 
first proton must be larger than the half width of the emis¬ 
sion of the second one”. However, the qualitative prediction 



FIG. 1. Schemes of different two-proton decays illustrating 
(a) true two-proton decay, (b) sequential decay, and (c) si¬ 
multaneous emission of two protons, but with two-body ^^F 
-I- p channel dominating. 


More detailed insight into the different decay modes can 
be found in the paper where the applicability of simple 
models of two-proton decay is carefully discussed. What is 
important here is that the reaction studied in the experiment 
of 0 is neither exactly within the scope of the original idea of 
Goldansky 0 for two-proton radioactivity nor exactly within 
the scope of idea of democratic decay [0-^ . This implies that 
other breakup mechanisms could also be important, and we 
now review such possibilities. 

Qualitatively different cases of two-proton emission are il¬ 
lustrated in Fig. 1^. Presumably democratic two-proton decay 
of the ^®Ne ground state (Fig. Q-^), where there are no 
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strongly correlated subsystems, is compared with decays of 
excited states in ^'*0 (Fig. 0) and i«Ne (Fig. |c, §), 
where binary decay channels are opened. The conceptual dif¬ 
ference between the experiments Refs. [0 and Q is that there 
is a sequential decay branch through the 1/2^ state at 2.37 
MeV in (Fig . |^b), which was shown to be important de¬ 
cay mode in ||m. In the case of ^®Ne sequential decay is not 
possible, but one-proton decay is still allowed to the weakly 
bound states in ^^F (5/2^ ground state, Eb = 0.6 MeV and 
1/2^, Eb = 0.11 MeV, Fig. ^). In the recent experiment 0 
the flux in the ^^F-|-p channel is orders of magnitude larger 
than in the three-body ^^0+p+p channel and can not be ne¬ 
glected in the study of the 2p channel. In such a case it will 
be very difficult to disentangle the three-body decay of '^^Ne 
states from the breakup of ^^F on a proton target (which also 
provides two protons in the final state), making the interpre¬ 
tation of the reaction more complicated. 

We therefore examine the nature of the reaction § by 
means of i?-matrix theory, to see whether the observed two- 
proton events can be explained in terms of three-body decay 
from the !“,£!*= 6.15 MeV resonance state in ^®Ne. We cal¬ 
culate the two-proton production cross section and branching 
ratios. 

Estimates. To estimate the cross section value for two- 
proton production via the resonance compound state with 
definite we use the standard formula (Ref. 0 ) 

(p\ — JL _ b'gFg _ 2J +1 _ 

{Er - Ey + FV4 (2Jm -b 1)(2J2„ -b 1) ’ ^ ’ 

where a and /3 label the entrance and exit channels, Jia and 
J 2 a are the spins of the particles in the entrance channel, F 
is the total (experimental) width of the resonance and Fi are 
partial widths. For the 1~ state F = 50 keV [|^. To obtain 
the complete cross section in the case of elastic scattering 
(a = /3) the Coulomb and potential scatterin^ogether with 
any interference terms must be added to Eq. (g). 

The resonant 2p production cross section via a 1“ state is 

<^2p [Er) j^2 p p ’ 

where, evaluated on the resonance, irlk^ = 0.31 barn, and 
rip(gs) is the decay width to the ^^F ground state. It can be 
seen from that cross section value of 310 mb is the upper 
limit for resonance processes via the 1“ state of ^®Ne at 6.15 
MeV. 

If the decay of the 1“ state to the excited 1/2’'" state in 
'^^F is negligible, the ratio rip(gs)/r is very close to unity. 
A more realistic estimate takes into account the variation of 
penetrabilities with decay energy. The conventional R-matrix 
formula is [0 

Flp(E) = 2 S(p E [E, ZC^F )), (4) 

where M/ is the reduced mass for particles with mass num¬ 
bers j and k. In this article the channel radius is varied 
from 2.5 to 4.5 fm to give an idea of the theoretical uncer¬ 
tainties associated with this parameter. Neglecting the dif¬ 
ference of spectroscopic factors Sip for the '^^F(gs)-bp and 
^^F(l/2’'')-bp channels we obtain rip(gs)/rip(l/2+) -- 2, and 
hence Fip(gs)/F ~ 2/3. If the inelastic channel '^^F(l/2’^)-bp 


has a large spectroscopic factor, it can influence drastically the 
cross section of the two-proton decay. This is connected with 
the fact that the large spectroscopic factor S^p will decrease 
the ratio Fip(gs)/r in Eq. (^) and will lead to smaller 2p pro¬ 
duction through the 1“ resonance. Unfortunately there was 
no experimental identification of the inelastic ^^F(l/2^)-bp 
channel in Q. Such identification is highly desirable to pro¬ 
vide more reliable estimates for the process. 

To estimate the simultaneous two-proton emission width 
we use a formula from ^ (similar to one given in 0): 


r2p(i52p) = 2 S2p 


nrUMl,r/^El^ 


- f' 


dx 


xPi, {xE2p,r., ZC^O)) Pi, ((1 - x)E2p, r^, zCo)) , (5) 


where h and h are the angular momenta of the 2 protons, 
and E 2 p is the three-body energy relative to the two-proton 
breakup threshold. The coefflcient in front of the penetra¬ 
bilities plays the role of a reduced width and is normal¬ 
ized in the spirit of the Wigner limit: without any barriers 
the width is the average inverse flight time for the distance 
rc/3. We certainly expect the two-proton spectroscopic fac¬ 
tor to be smaller than the one-proton spectroscopic factor: 
S 2 p < Sip. However, assuming S 2 p = Sip the upper limit 
for the estimated branching ratio with H = 0 and Z 2 = 1 is 
r 2 p/r = (1.5 — 5) X 10“"' (for rc in the range 2.5-4.5 fm, see 
solid curve in Fig. ^). 



FIG. 2. The ratio F 2 p/r as function of channel radius 
Tc estimates in two-proton and ‘diproton’ models. Solid 
and dashed curves correspond to 1“ state (simultaneous and 
‘diproton’ emission) and dotted curve corresponds to 2“ state 
(simultaneous emission only). 

Another simple estimate of the branching ratio is provided 
by the diproton model: 

r2p{E2p) = 2S2p Pi{E2p - e,r,,2ZC<^0)) , (6) 

where I = 1 for the 1“ state and the average internal energy 
of the ‘diproton’ e ~ 0.3 MeV is estimated from the experi¬ 
mental distribution Fig. 4 in The usual assumption is to 
take e = 0 (see for example [ iflUT^ ) or e = 0.05 — 0.1 MeV 
(for example 0J0) as the energy of the ‘virtual state’ in 
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two-proton system, but these prescriptions run into method¬ 
ological problems [EOIH . There is no need to follow these pre¬ 
scriptions if one has an idea what this energy actually is. The 
use of the ‘experimental’ diproton energy in Eq. (^) provides 
an estimated branching ratio r 2 p/r = (2 — 10) x 10““* which 
is in a very good agreement with the simultaneous emission 
estimate above. Figure ^shows the ratio r 2 p/r for both types 
of estimates as a function of channel radius r^. 

Discussion. In the estimates above we have made several 
assumptions. Each of them is likely to give us the upper limit 
for the 2p production cross section. We should also note that 
the source of the theoretical uncertainty in the 2p cross sec¬ 
tion is a variation of the channel radius re- The source of 
experimental uncertainty in the two-proton production cross 
section is connected with different assumptions about the de¬ 
cay process (‘diproton’ decay or ‘democratic’ decay) in the 
analysis of the data, because of the limited acceptance of the 
experiment. 

The estimates for r 2 p/r ratio are consistent with theoreti¬ 
cal estimates in 0. We deduce from Eq. (|^ an estimated two- 
proton production cross section through the 1~ resonance of 
0.03 — 0.2 mb, much smallei]^ than the measured value 1.5 —4.0 
mb reported in |^. 



FIG. 3. Schematic presentation of (a) emission of two pro¬ 
tons via the decay of resonance in ^®Ne, and (b) dominating 
mechanism of breakup on protons. 

There are two possible explanations of the large two-proton 
production cross section seen in the experiment of [Q|. 

(i) According to the level scheme in the isobaric nu¬ 
cleus, and to the experimental evidence from the 2~ 

and 3“ states are located only slightly higher than the 1“ 
state. Two-proton contributions from the 3“ state should be 
negligible, as = 1 and I 2 = 2 for the simultaneous emis¬ 
sion model Eq. (^, or i = 3 for ‘diproton’ emission Eq. (|^. 
The ‘diproton’ emission from a 2~ state is parity forbidden as 
has been mentioned in [|^ , but simultaneous two-proton emis¬ 
sion is allowed with the same quantum numbers = 0 and 
^2 = 1 as from the 1“ state, with a width comparable with the 


^ The branching ratio only gives the experimental cross sec¬ 
tion if it multiplies a cross section of 3660 mb |^, in note 16] 
for Ip resonance production from Eq. (|^. 


‘diproton’ mechanism (see Fig. |^, dashed line). Using Eqs. (P) 
and (^ we obtain the value of ~ 0.1 — 0.4 mb. This value 
would decrease significantly the difference between measured 
and estimated two-proton production cross sections. 




FIG. 4. Qualitative comparison of different possible contri¬ 
butions to the two-proton events, (a) Shows the estimated 
contributions of possible breakup events in experiments at 
33 and 44 MeV beam energies (for different channel radii). 
Fig. (b) compares the expected shape of the strength func¬ 
tions for breakup events (solid curve is the same as in up¬ 
per panel), resonance events (from 1“ and 2~ states; dotted 
curve), and the experimental distribution approximated by a 
gaussian (dashed curve). 

(ii) The other possibility is that most of the two-proton 
events are actually coming from the breakup of ^^F on pro¬ 
tons, rather than decay of resonances in ^®Ne. Figure ^ 
schematically outlines the difference between the two pro¬ 
cesses. It is difficult to estimate the breakup process con¬ 
sistently, as it can involve the complicated interplay of two- 
body and three-body dynamics. What is important here is 
that this contribution is not limited by the form of Eq. (^) 
which applies only for 2p production via a resonance. 

The energy dependence of the two-proton breakup chan¬ 
nel can be roughly estimated by the formula for simultaneous 
two-proton emission, Eq. (^, as both protons have to pen¬ 
etrate through the Coulomb barrier. However, the channel 
radius no longer has a well defined meaning, although it is 
certainly expected to be larger than the values used in esti¬ 
mates for the decay of a compound state. Imagine for the 
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moment that there is no resonance contribution to the two- 
proton production at all. The estimated energy dependence 
of the breakup cross section is shown in Fig. ^a. In the type of 
experiment described in [Q| , all beam energies below maximal 
are present with comparable probabilities because of the thick 
target used. The estimated ratio of breakup events observed 
in experiments with 44 MeV and with 33 MeV beams is 20- 
40 depending on channel radius. This is given by the ratio of 
integrals of the intensity shown in Fig. ^ up to Er = 1.84 
MeV and Er = 1.23 MeV respectively. The main contribu¬ 
tion of two-proton events from inelastic breakup comes from 
the maximal energies available in the experiment. In Fig. 
^ this contribution is qualitatively compared with the ex¬ 
pected contribution from 1“, 2“ states and the distribution 
of actually observed events, which is broad due to the energy 
resolution being low in the experiment. In our view all the 
above mechanisms could contribute to the experimental cross 
sections. 

The ratio of the two-proton cross section for the 
EiabC^F) — 33 MeV measurement compared to the 

EiabC^F) = 4:4 MeV measurement could be possible evidence 
for identification of decays via the 1“ resonance, and the ex¬ 
periment 1^ sees a 7-10 fold suppression at the lower en¬ 
ergy. However, we have to be careful that the growth of direct 
breakup does not dominate any resonance contributions. No 
resonant 2p contribution is expected in the 33 MeV measure¬ 
ment. The estimates of Fig. W shows that if the number of 
two-proton events in 33 MeV measurement is only 10 times 
lower than in 44 MeV measurement, then extrapolation of 
the intensity of breakup to higher beam energy gives more 
than enough 2p events to explain the whole 2p intensity in 44 
MeV measurement. This indicates that further experimental 
evidence will be need to discriminate between direct breakup 
and resonant decay. 

Conclusion. We have shown that it is likely that the two- 
proton decay of the 1“ state is sufficient to explain only a 
small fraction of the 2p events reported in Ref. Q. It is plau¬ 
sible that the balance of 2p events is connected either with 

(i) excitation of a 2~ state located a little higher than the 1“ 
or/and with (ii) the nonresonant breakup of ^^F on protons. 
Identihcation of the inelastic breakup channel ^^F(l/2"'') +p 
for the states involved in the two-proton emission is also de¬ 
sirable for a refined interpretation of the data. A complete 
kinematics experiment would be the most useful in the study 
of 2p emission from '^^Ne. It would allow the energy of the 
decaying states in ^®Ne to be fixed, and would thus make 
interpretation of the the experiment clearer. 
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